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GEOMETRIC  RETRIEVAL  PROBLEMS 

Richard  Cole  and  Chee  K.  Yap 

Courant  Institute  of  Mathematical  Sciences 
New  York  University 

Abstract 


A  large  class  of  geometric  retrieval  problems  has  the  foUowing  form.  Given  a  set  X  of 
geometric  objects,  preprocess  to  obtain  a  data  structure  D{X).  Now  use  D(X)  to  rapidly  answer 
queries  on  X.  We  say  an  algorithm  for  such  a  problem  has  (worst-case)  space-time  complexity 
0(f(n),g(n))  if  the  space  requirement  for  D(X)  is  0(f)  and  the  'locate  run-time'  required  for  each 
retrieval  is  O(^).  We  show  three  techniques  which  can  consistently  be  exploited  in  solving  such 
problems.  For  instance,  using  our  techniques,  we  obtain  an  0(n~'^',lognlog(l/£))  space-time  algo- 
rithm for  the  polygon  retrieval  problem,  for  arbitrarily  small  t,  in^oving  on  the  previous  solution 
having  complexity  0{n',\ogn). 
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1.  Introduction 

We  consider  the  class  of  geometric  retrieval  problems  described  in  the  abstract,  above.    The 
following  are  typical  examples,  many  of  which  have  been  extensively  studied. 

1.  Point  Location  [LT,P1,P2,LP,K]:  X  is  a  set  of  polygonal  regions  which  form  a  partition  of 
the  plane;  a  query  is  a  point  x,  and  the  answer  is  the  region  containing  x. 

2.  Nearest  Neighbors:  X  is  a  set  of  points;  a  query  is  a  new  point  x,  and  the  required  answer  is 
the  set  of  points  which  are  nearest  to  x. 

3.  fixed  Radius  Nearest  Neighbors  [BM,P2]:  Let  r>0  be  a  fixed  parameter  of  the  problem.  X 
is  a  set  of  points,  a  query  is  a  point  x,  and  the  answer  is  the  set  of  points  in  X  at  distance  ^r 
fromx. 

4.  Polygon  Intersection  [EKM]:  X  is  a  set  of  polygons;  a  query  is  a  polygon  P,  and  the  answer 
is  the  set  consisting  of  those  polygons  whidi  intersect  P. 

5.  Polygon  Retrieval  [EKM,W]:  Let  A:  be  fixed.  X  is  a  set  of  points;  a  query  is  a  it-gon  P,  the 
answer  is  the  subset  of  X  contained  inside  P.    This  is  a  generalization  of  the  well-known 
'range-query'  problems  where  the  query  is  a  rectangle  (the  'range')[FB,BS,Lu]. 
Generally,  the  efficiency  of  the  solution  is  measured  by  three  parameters:  the  space  taken  to 

store  the  preprocessed  structure  D{X),  the  time  taken  to  answer  a  query,  and  the  preprocessing 
time  taken  to  construct  D{X).  We  omit  reference  to  the  third  of  these  mainly  for  the  sake  of  sim- 
plicity. Also  in  some  applications  the  preprocessing  time  (a  one-time  cost)  is  less  important  than 
the  retrieval  time  or  the  space  to  store  the  preprocessed  structiu-e  (eg.  flight  simulation  graphics). 
It  is  usually  relatively  easy  to  determine  the  preprocessing  time  for  our  algoritlims;  typically,  it  is 
0{n.g{n)),  where  g{n)  is  the  space  complexity.  We  use  the  'locate  run  time'  in  the  sense  of  [W] 
for  our  qtiery  time:  it  is  the  time  to  locate  the  first  object  in  the  answer,  where  it  is  assumed  each 
additional  object  can  be  output  in  constant  time. 


Our  initial  motivation  was  the  triangle  retr:e\'al  problem.  It  can  be  used  to  solve  the  polygon 
retrieval  problem,  by  decomposing  the  polygon  into  disjomt  triangles.^  [EBOil]  obtained  a  space- 
time  bound  of  0(n' ,logn)  for  this  problem;  they  ware  cpnsidepug  the  more  general  problem  of 
polygon  intersection  but  the  space  complexity  of  triangle  retrieval  jdominatr:  in  their  solution. 
They  posed  the  challenge  of  reducing  this  space  con^lexity. ,  We  stiow  that  the  space-time  com- 
plexity can  be  drastically  reduced  to  (Xn-^^'.lognlogCl/c))  for  any  e>0;  we  also  obtain  a  solution 
with  complexity  0(n-,lognloglogn).  Both  solutions  admit  extensions  to  higher  dimensions. 

An  alternative  solution  to  the  triangle  r^trie^!^  probl^  was  giypi  in  [W]  It  has  complexity 
0(n,n''^.  This  solution  has  since  been  iiaproysd^j-Edel^hpijgier^tQLa-compl'^ri^^  of  0{n,nP^) 
(private  communication  from  Willard).    Tt  h^s  leeeaitly^ij^e^.ert^^  to  H-mension  3  by  Yao 

[Yl],  but  the  method  apparently  breaks  down  for  dimensions  ^5.  j_  pj.-, 

I 
Our  algorithms  are  based  on  three  tediniqi^K.  They  are  (|uite  robust  in  that  they  generalize 

to  higher  dimensions  and  can  be  applied  to  a  variety  of^rdat^  j^oblems.  We  illustrate  these  tech- 
niques  on  the  triangle  retrieval  problem.  The  first  tedmique  is  to  potion  the  query  space  (if 
there  are  r  parameters  for  the  query  we  can  think  of  thi^  as  an  r-dimensional  space)  into  regions 
for  which  some  property  of  the  answer  is  invariant.  Tlis  simplest  example  of  such  a  partition  is 
where  each  region  has  a  constant  answer.  This  technique  is  known  by  various  names  in  the  litera- 
ture, eg.  the  'locus  method'  [BM],  the  'priority  approach'  [SSS,Y2].  The  second  technique  is  to 
take  advantage,  where  applicable,  of  the  fact  that  invariant  properties  of  the  answers  are  similar  in 
nearby  regions,  and  storing  the  shared  information  only  once.  In  computational  graphics  this  is 
known  as  'coherence'  [SSS].  The  third  technique  is  a  form  of  divide  and  conquer:  we  decompose  a 
large  difficult  problem  into  one  small  difficult  problem  and  several  large  easy  problems. 

In  the  next  three  sections  we  illustrate  these  techniques  by  obtaining  a  solution  of  the  trian- 
gle retrieval  problem.  In  the  fifth  section  we  provide  a  solution  to  the  nearest  neighbors  problem: 
this  is  the  fixed  radius  nearest  neighbors  problem  except  that  now  the  radius  is  also  a  parameter  of 
the  problem.   In  the  sixth  section  we  give  a  solution  for  tetrahedral  retrieval  (in  3  dimensions).  In 


the  seventh  section  we  briefly  discuss  other  problems  to  \Vhi(S  we  apply  these  techniques.  In  the 
final  section  we  provick  a  summary  of  GtoJresiilter-' On-thdi^y'tO  a  solution  of  the  triangle  prob- 
lem we  solve  the  half-plane  problem:  "vtlaSi  poinlf-l'e-lb'the^eft  (right)  of  a  given  line?',  and  the 
wedge  problem:  'which  points  lie  betWeein  two  temi-infinite""tifies  Z-j  and  Lj.  whidi  meet  at  their 
common  endpoint?';  in  both  problems  the  pcSiffe  lieon  ek  ^lane. 

2.  Half  Plane  Retrieval  -'  ,-:  lb  ;  a  -i-         '    "= 

We  illustrate  the  first  technique  by  showiiig  how 'iR>  solve  Afr  hMf-plane  problem.  We  draw 
a  vertical  baseline  to  the  right  of  the  ?)iMateHitA7-€©ad<ief  dief^qu^'ffife  iL;  it  cuts  the  baseline  in 
some  point  P.   Comider  rotating  t 'about- P^3^^i*ffl^il*feftied  alt Ihfe  ^©iite  Associ- 

ate with  P  a  circular  ordered  list  of  points  in  X.       •'--'  ^nciansn-;;  ui  r-  .. 

We  can  form  mtervals  on  the  baseline  consisting  of  those  points  P  whose  associated  lists  are 
identical.  The  separators  of  these  intervals  are  the  intersections  of  the  baseline  with  lines  through 
pairs  of  points.  Thus  there  are  O(n-)  intervals,  and  with  each  interval  we  associate  a  list  of  n 
points. 

To  determine  which  points  lie  to  the  left  (right)  of  L,  we  find  in  which  interval  P  lies,  and 
then  search  the  corresponding  list  (stored  in  an  array)  to  find  those  points  to  the  left  (right)  of  L. 

In  fact,  this  solution  solves  the  based-wedge  problem  also:  Vhich  points  lie  between  lines  LI 
and  LI,  where  L\  and  LI  meet  on  the  baseline?'. 

This  solution  has  complexity  0(n^,log/:).  Below  we  give  a  version  of  our  solution  having 
complexity  0(n-,logn). 


For  the  half-plane  problem,  Chazelle  (private  comrmipication)  and  Gtaibas  (private  communi- 
cation) have  independently  shown  that  'duality'  can  be  exfrfoited  to  give  a  bound  of  O(n,log/i). 
However,  their  method  does  not  seem  to  prov^Je  an-efficiej;it  splutioaof  the  based- wedge  prob- 
lem. ,,  ,^ ,_  V  : 

3.  The  Half  Plane  Problem  Again 

We  illustrate  the  second  technique  on  the  half-plane  piroblem.  Observe  that  the  lists  of 
points  cxjrresponding  to  adjacent  intervals  differ  on  exactly  two  points;  so  it  seems  excessive  to 
store  the  whole  list  twice.  /  I  ;       ~ 

In  [Co]  the  following  problem  is  solved.!  Qven  a  sequencejaf-iirtfl'dered  lists  of  n  items 

.  ■        I 

each,  such  that  the  sum  of  the  differences  between  adjacent  lists  is  A,  we  can  store  these  lists  in 

space  0{n  +  h),  while  allowing  the  followjig  queries  to  be  carried  out: 

a)  report  the  ith  element  in  the  jth  list  in  tim^O,flog('l-bW)     -  ■      _,    , 

b)  report  the  items  between  the  ith  ar'Ajih  m  t|ielA-li8t  in.tiij»^(Iog(«-i-.'!)+j-i). 

It  is  often  possible  to  carry  out  a  variant  of  opera^ij^i  QJ),  where  ihc  values  of  i  and  j  are  not 
given  explicitly,  as  is  the  case  for  this  application.  We  apply  this  solution  to  our  problem  of  stor- 
ing O(n-)  lists  of  length  n,  where  adjacent  lists  differ  on  exactly  2  points.  We  obtain  a  solution  of 
complexity  0{n-,\ogn). 

4.  The  Wedge  and  The  Triangle 

We  illustrate  the  third  technique  first  by  showing  how  to  solve  the  wedge  problem  in  the 
plane,  and  second  by  showing  how  to  solve  the  triangle  problem.  WLOG  we  assume  that  wedges 
are  convex. 

A  first  attempt  at  the  wedge  problem  would  be  to  proceed  as  in  section  2,  except  now  we 

f(xa.  a  list  for  each  point  in  the  plane.    Regions  of  points  with  the  same  list  are  found.    These 

regions  are  defint.  1  by  the  lines  through  pairs  of  points;  there  are  0{n*)  of  them.  It  leads  to  an 


0(n^,logn)  solution.  (It  can  be  improved  usii^  the  methed  of  section  3  to  0(/i'',logn)).  Instead, 
divide  the  set  of  n  points  by  some  vertical  line  Ly  into  two  s6ts  oi-n/2  points.  Let  L^,  L^  be  verti- 
cal lines  to  the  left  and  right  of  all  the  poiitts,  attd  c^  thft-TegibiH  between  these  lines  R^,  /?^. 
WLCXj  suppose  the  comer  of  wedge  W  is  in  /?^.  If  W  lies  completely  within  R^  it  is  a  wedge 
query  on  n/2  points;  otherwise  there  are  two  possibilities,  illustrated  below. 


■>i-:si'-.3c  r  n-. 


Ziii.  Jnoaqbf    .35W?r  ^ti 


Ether  (a)  the  result  of  W  is  the  union  of  the  results  of  three  queries:  VV3,  a  wedge  query  on 
Xr\Ri,  W,  and  W,,  based-wedge  queries  on  XQRjj,  wmTw^f  basfetf  On  Lr  arid  W,  based  on  L^^.  Or 
(b)  the  result  of  W  is  the  union  of  two  qu^ri^  W^,^^  i»ix^  ^6ry  on  xn/?^;  i/,  a  half-space 
query  on  Xn/?^.  Using  recursion,  and  by  the  results  of  secttOM  2  and  3,  this  yields  a  solution  of 
complexity  C)(/i-,log-n). 

Instead,  suppose  we  split  the  set  of  points  into  n'  sets  of  n^~*  points,  and  these,  in  turn  into 


7^-    sets    of    n'   ^'^ 


and    these    in    tiun    into    n 


a'2/. 


sets    of    size    n 


where 


Jk,  =  t(H-3/2-(-(3/2)--H...-i-(3/2)'),  for  fc,<l,  ...  .  This  yields  a  solution  having  complexity 
0(n-*«,log/ilog(l/£)). 

Alternatively,  we  could  start  by  dividing  the  set  of  points  into  2  sets  of  «/2  points,  then  ..., 
each  set  of  nid  points  into  d^'^  sets  of  n/(d)-'''^  points,...  .  This  yields  a  solution  having  complexity 
(9(n^,lognloglogn). 

In  computing  the  complexity  in  these  latter  cases  one  should  note  that  if  a  set  K  is  divided 
into  equal  sized  disjoint  sets  Y.,Y-,,...,Y^,  we  have  to  be  able  to  answer  based  wedge  queries  on  all 


the  sets  \jY,and\j  Y„  l<i<Jk.  TTiat  is  0(k)  sets  of  si?e  Oi\Y^,. 


s'\  s-i 


The  triangle  query  can  be  decomposed  in  a  similar  way: 


giving  one  wedge  query  W  and  a  a  2-cor/Mt  qitoy^CVora  tifemglequCTy  on  a  point  set  of  half 
the  size.  The  2-comer  query  can  be  decomposed-in  a  siiiiifiarwayi'tilsb:'^' 


it  is  either  the  union  of  2  based  wedge  queries  W-  and  W-,  plvs  a  2-comeT  query  C  on  a  point  set 
of  half  the  size,  or  it  is  the  union  of  a  wedge  query  W  plus  a  1- comer  query  C.  The  1 -comer 
query  can  be  decomposed,  also: 


it  is  the  union  of  2  based  wedge  queries  W-^  and  W,  plus  a  1-comer  query  C  on  a  point  set  of  half 
the  size.   On  analysing  this  we  find  we  have  a  solution  of  complexity  C>(n',log'n). 

Again,  we  could  divide  the  n  points  into  n*  sets  of  n^~*  points  each,  etc.   We  then  obtain  a 
solution  with  complexity  0(n'*',lognlog(l/£)).   Or  we  could  divide  sets  of  size  n/d  into  d^'^  equal 


sized  sets,  yielding  a  solution  with  complexity  O(n',lognloglogn); 

:\    A'    jii 
5.  Disk  Retrieval 

The  query  is:  Svhich  points  lie  inside  the  disk  with  radius  rjxntaed  at  a  given  point?'.  This 
is  a  generalisation  of  the  fixed  radius  nearest  neighbor  problem  [BM], 'where  the  radius  r  is  fixed 
beforehand.  It  also  solves  the  k-ntaie&t  neighbor  problem  [L],  improving  on  the  0(nk)  space 
complexity  of  [L],  for  ;k  =  co(log'n)  (the  problem  is  to  list  the  t  closest  neighbors  of  a  given  point). 

We  give  a  method  having  complexity^|^lQ^p,l(^5i),,  It-jhould.be^observed  that  the  count- 
ing of  the  running  time  differs  from  the  other,  gages.,,  ^^i^^we  ^orti^.the  search  time  over  the 
reporting  time;  that  is,  if  the  reply  cosntains  s  pointy,  the  search  time  is  0(s  +  \ogn).  This  is  the 
same  counting  method  as  is  used  by  Qiazelle'lCWinlhis  paper  on  filtered  search. 

We  use  a  recursive  data  structure,  whidLwe  think  of  as  a  recursive  tree.  A  node  in  the  tree 
handles  queries  on  some  subset  S  of  the^pointSi  IrHaTa  local  data  structure  able  to  handle  queries 
having  fewer  than  log-n  points  in  5.  Also,  it  has  a  local  data  structure  to  enable  one  to  determine 
in  O(logn)  time  whether  there  are  at  least  \og-n  points  tfiiiS  5  iil  the  answer.  The  local  data  struc- 
tiu'es  are  described  later. 

The  root  of  the  tree  handles  queries  on  aU  n  points.  If  a  node  of  the  tree  handles  queries  on 
a  set  5,  with  |5|3:log-n,  then  it  has  two  diildren  which  handle  queries  on  the  disjoint,  equal-sized 
sets  5;,  5;  with  5- U52  =  5.  We  show  below  that  the  local  data  structirre  for  a  node  handling  a  set 
of  m  points  requires  space  O(mlog'n).   Thus,  the  overall  space  requirement  is  O(nlog'n). 

Next,  we  describe  the  search  algorithm.  Suppose  we  are  searching  a  node  N  of  the  tree, 
handling  queries  on  a  set  S  of  points.  K  there  are  at  least  log-n  points  from  5  in  the  answer,  then 
we  search  both  children  of  N.  Otherwise,  we  use  the  local  data  structures  at  N  to  list  those  points 
from  5  in  the  answer. 

We  establish  the  running  time.  We  show  below  that  for  a  local  data  structure  handling  a  set 
S,  in  time  0(log/i)  we  determine  whether  there  are  at  least  log-n  points  from  5  in  the  answer,  and 


in  time  0(logn  +  j)  we  list  the  s  points  from  S  in  the  answer,  for  s<\og-n.  We  charge  the  0(s) 
work  at  a  rate  of  0(1)  to  eadi  point  listed.  The:<?OQg^)  wofk  is  charged  to  the  parent  node, 
from  which  the  recursive  call  was  made.  The  charge;  aJ  ea^  qpde  iS; evenly  divided  between  the 
log-n  nearest  neighbours  among  the  poipts  hwdJjBdcat  that  node.-rThis  is  a  charge  of  0(lAogn) 
per  point,  every  one  of  which  lies  in  the  ejisw^y:  jEach  point  receives  sucfe  a  charge  at  most  once 
per  level  of  the  recursive  tree.  There  are  Qtdy  €f(}ogfi)  levels  in  the  tree,  and  so  there  is  a  charge 
of  only  0(1)  per  point  in  the  answer.  There  is  also  a  change  of  0(logn)  for  Jie  search  at  the  root. 
So  the  overall  running  time  for  a  seardi  'S  0(\ogn^i^i  whore  s  is  the  number  of  points  listed. 

Now,  we  discuss  the  local  data  structures.  We  build  a  log^n-nearest  neighbor  diagram  -  each 
region  in  this  diagram  has  the  same  lo^-n . 'A iKarest  neighbor..  AiUttle  thought  fibows  that  the  k- 
nearest  neighbor  diagram  is  a  coarsening  of  the.Sttpe?pQsto»dnnQfi.^  (^T.l)".Voronoi  diagram  and 
the  Jt-Voronoi  diagram  (see  [L]  for  definitf—.naad  sonStnidiOTO-ii'n»  fe\*Wonci  diagram  has  size 
Oink).  On  examining  the  method  for  constnicting^^hRirfc-Vptonoi  diagram  from  the  ()t-l)- 
Voronoi  diagram  we  notice  that  the  z>i^^  of  the  two  tdiagramS' oidyristerserf:  at  their  endpoints 
[L];  thus  the  size  of  the  ik-nearest  neighbrji  diagrcjnaSjfibC^fc)*  :  r  >     .— 

We  use  the  log-n -nearest  neighbor  diagram  tw  determine  whether  there  are  fewer  than  log-n 
points  in  the  answer.  To  handle  those  queries  with  fewer  than  log-n  points  in  the  reply  we  also 
build  the  ifc-Voronoi  diagram,  for  fc  =  logn,21ogn,41ogn,..., log-n,  storing  a  list  of  the  it  nearest 
points  for  each  region  in  each  diagram.  To  respond  to  such  a  query  we  start  by  searching  the 
logn-Voronoi  diagram.  Having  found  the  correct  region  for  the  query,  we  then  search  the  list  of 
points  associated  with  that  region.  If  some  of  the  points  do  not  lie  within  the  query  our  reply  con- 
sists just  of  those  points  that  do  lie  within  the  query,  among  the  points  in  this  list.  Otherwise,  we 
proceed  to  search  the  21ogn-Voronoi  diagram  in  the  same  way,  and  then  the  41ogn-Voronoi 
diagram,  and  so  on,  until  we  find  a  point  that  does  not  lie  within  the  query.  Any  such  search 
requires  time  O(logn-^j),  where  s  is  the  number  of  points  in  the  reply. 


We  need  a  list  of  length  k  for  each  region  in  the  )t-Voronoi  diagrams  that  we  store.  This 
leads  to  a  space  requirement  of  0{n\og*n)i  We  sHOW  hOW  to  reduce  this  to  0(/ilog-n),  without 
affecting  the  order  of  the  search  time.  Notice  that  the-lfets^feOTi«sponding  to  adjacent  regions 
differ  on  exactly  one  point.  So  if  we  coilld  fonn  a  sequence  oP-is  te^m,  with  neighbors  in  the 
sequence  being  adjacent,  then  these  k  regions  would  have  at  most  2Jtdistiiiict  points  among  their  k 
nearest  neighbors.  We  would  then  store  in'<»le  list  the  ap  to  2*  points  fw  these  k  regions.  When 
searching  the  list  we  would  now  alsokftep  a  count  of  h©w  many  points  were  found  inside  the 
query;  this  would  enable  us  to  inroceed^^^W  above.  -R  'remtos  to  show  how  to  form  these 
sequences.  ' -Isn  3d- 7r-.r'soU  bliiid  3^7    .^a:^.-, 

We  form  the  dual  of  the  region  graph.  Tliit'is,  eadlfyegioa-iSLwpIaoKl  by  a  vertex,  and  ver- 
tices are  joined  by  an  edge  if  and  only  if  the  coftttepondbiiafBgiojB  are  adjacent.  We  build  a  span- 
ning tree  of  the  dual  graph,  and  then  tatoe'an  EGTeriDtfcaiirotihe  8panr.'.ig  tree.  If  the  graph  had 
r  nodes  then  the  length  of  the  circuit  is  £>(/■) ,'«id-iiirfieirt6verj(nnode- occurs  on  the  circuit.  On 
replacing  the  nodes  by  the  corresponding'- iSgidas  We  lQ|Kt)ia 'se(^ifi»ice  o£length  0(nk)  for  the  k- 
Voronoi  diagram.  We  divide  this  sequence  into  0(n)  sfequ^ces,'  each' of  length  k,  giving  us  the 
desired  sequences  of  adjacent  regions.  The  overall  storage  requirement  is  reduced  to  0(nlog^«). 
This  completes  the  data  structure  for  handling  the  queries  with  small  answers. 

As  claimed,  it  takes  O(logn)  time  to  determine  whether  there  are  at  least  log-n  points  from 
5  in  the  answer.  Also,  for  a  node  from  which  no  recursive  call  is  made,  it  takes  time  0(s  +  \ogn) 
to  list  those  s  points  from  5  in  the  answer. 

The  above  method  works  regardless  of  how  the  n  points  are  divided  into  two  sets  of  n/2 
points.  We  show  this  division  can  be  chosen  to  our  advantage,  and  thereby  reduce  the  space 
requirement  to  O(nlog-n),  without  affecting  the  running  time. 

For  each  point  set  of  size  m  obtained  by  the  recursive  dividing  we  keep  a  clogm-nearest 
neighbor  diagram,  for  some  constant  c,  rather  than  a  log-^n -nearest  neighbor  diagram,  as  before. 
To  provide  a  search  time  of  0{]ogn  +  s)  we  require  that  each  set  of  clogm  nearest  neighbors  be 
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(roughly)  evenly  divided  between  the  two  sets  of  m/2  points. 

Let  5  be  the  set  of  points  handled  by  an  intern^- tfodi^V  ©f  the  tree,  and  let  X  be  the  set  of 
points  handled  by  its  parent.  For  N,  with  respect  to  a  location  /,  we  say  the  new  points  are  those 
points  among  the  clogm  nearest  neighboiifs -irf- J/  in-S  that  are  not  fflnong  the  clog (2m)  nearest 
neighbours  of  I  inX  {\X\  =  2m).  Tliere  are  0(lffgm)  tiew  points.     " 

We  charge  the  work  as  before,  with  one  difference.  The  work  charge  to  a  node  is  divided 
among  the  new  points,  which  is  an  0(1)  charge  per  new  point,  each  one  of  which  lies  within  the 
answer.  With  respect  to  a  given  location,  eadi  poiiit  is  a  new  pioint  at  Kiactly  one  node  of  the 
tree.   So  again,  we  obtain  a  running  time  of  O^o^rt'^sy,  >Vh^  jls'thetluftdser  of  points  listed. 

We  show  such  a  division  of  a  set  of  n'pbints  info  t^o  sets  6Pn/2  pbints  exists;  in  fact,  we 
show  that  with  high  probability  a  random  ^vision  is 'su^"a(i^da'.  tb'dd  so  we  solve  the  fol- 
lowing more  general  problem.  Given  sets  S,  bf'aze'  cioin^^i=lX\'.':,dniogn,  all  chosen  from 
some  set  5  of  n  vertices,  for  some  constant  d,  vmire'd  depeads  one.  We  want  a  division  of  5 
into  two  disjoint  sets  A,  B,  of  nil  vertices  each,  such" that  l,J,(iAf  a  alogn,  |5,nfl|  a  alogn,  for 
j  =  l,2,...,^logn,  for  some  constant  a,  0<a<l. 

By  counting  the  number  of  choices  of  A  that  violate  one  of  the  conditions,  we  find  that  most 
choices  of  A  satisfy  the  conditions,  for  large  enough  c,  small  enough  a,  and  large  enough  n.  This 
proves  that  the  construction  can  be  made,  in  a  probabilistic  fashion  at  least.  We  do  not  include 
the  details  of  the  counting  argument  as  they  are  straightforward  and  tedious.  This  construction 
suggests  a  further  problem:  can  it  be  done  deterministically? 

6.  Tetrahedral  Retrieval 

The  query  is:  Svhat  points  lie  inside  a  tetrahedron?'.  The  points  lie  in  3  dimensional  space  in 
this  problem.  We  first  show  how  to  do  half-space  retrieval  in  3  dimensions,  then  extend  the  result 
to  two  special  types  of  retrieval:  based  cone  retrieval  and  aligned  tetrahedral  retrieval,  defined 
below.  Fmally,  we  use  these  to  solve  the  tetrahedral  retrieval  problem. 
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Half  Space  Retrieval 

The  half-space  retrieval  problem  consists  gf  th§  foljowing  query:  Svhat  points  lie  to  one  side 
(the  left  or  the  right)  of  a  given  plane?'.  As  with  thphftjf -plane ^yptrieval  problem  we  consider 
rotations.  For  this  problem  we  fix  a  base  plane  (z=0,  s^y)_  apd  then  consider  rotations  of  planes 
about  their  intersection  with  the  base  plane  (a  line).  Wft  seek  'regions'  of  lines  which  produce  the 
same  rotational  order  for  the  points. 

Consider  a  pair  of  points  p  and  q.  Let  the  intercept  Qf  the  line  Jhrough  p  and  q  with  the 
base  plane  be  the  point  r.  Let  L  be  the  line  in  the;,  base  plane  abpftt  which  we  are  rotating  the 
plane.  We  see  the  rotational  order  of  the.;pcane?j?vas(4'vdq5«id*:SDlely  on  v/hich  side  of  L  the 
point  r  lies.  So  this  becomes  almost  a.half-plan^  ^^trieyaj  problem  on  n^  points,  the  points  that 
are  the  intercepts  with  the  base  plane  of  lines  thj'ough  pairs,  of  ^the  original  points.  We  obtain  n* 
intervals  on  a  base  line  (v=z=0  say),  and  with  e^ch  interval  there.is  an  associated  list  of  n~  inter- 
cept  points.  To  each  initial  segment  of  every  sudi  list  there  corresponds  a  list  of  n  points,  which  is 
the  rotational  order  of  the  original  points  with  resp^  to  rotatipns  about  an  appropriate  line:  any 
line  that  separates  the  initial  segment  of  the  intercept  points  from  the  remainder. 

We  use  the  method  in  [Co]  to  store  the  lists  of  intercept  points,  storing  them  in  a  table  T  in 
space  0{n*).  Next,  we  observe  that  there  are  Oin'^)  sets  into  which  the  intercept  points  can  be 
divided  by  lines  in  the  base  plane,  and  that  every  entry  in  T  corresponds  to  a  unique  partition  of 
the  intercept  points.  Further,  for  every  such  partition,  there  is  an  entry  in  T.  It  remains  to  store 
a  list  of  length  n  with  each  entry  in  T,  namely  the  rotational  order  of  the  original  points 
corresponding  to  that  partition  of  the  intercept  points.  We  notice  that  lists  corresponding  to  'adja- 
cent' regions  (of  lines)  differ  on  exactly  two  points.  Further  the  sets  of  'line  regions'  form  a  con- 
nected graph.  So  by  taking  a  spanning  tree,  etc.,  as  described  above,  we  can  store  the  0{n*)  lists 
of  length  n  in  space  0{n*),  while  allowing  queries  to  be  answered  with  an  overhead  of  only 
0(log/i)  time.  The  overall  complexity  is  0{n*,\o%n).  We  are  not  interested  in  the  list  of  intercept 
points  per  se;  they  merely  provide  a  way  to  determine  which  list  of  the  original  points  to  search. 
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achieved. 

Based  Cone  Retrieval 

A  based  cone  query  is  the  intersection  of  three  half  spaces,  each  defined  by  a  plane,  where 

the  intersection  of  the  three  planes  (the  comer)  lies  on  a  given  base  plane.   We  need  to  consider 

another  problem,  the  three  dimensional  wedge  query.   This  is  the  intersection  of  two  half-spaces 

defined  by  planes.    We  solve  it  by  iterating  the  3-dimensional  half-space  retrieval  solution.   By 

considering  the  dual  problem  for  the  half-space  retrieval  problem  (points  are  mapped  to  half- 

spaces  and  half-spaces  to  points),  we  see  there  are  only  £?(rt^)  different  solutions  to  the  half-space 

problem  (each  solution  corresponds  to  a  region  in  the  dual  problem).  For  each  of  these  solutions 

r  ..'I    .J  ;,  J-  7C  L;r  r.-rJ-j;       -.^  .rr:. 
we  construct  a  data  structure  allowing  half-space  queries  to  be  performed.   This  provides  a  data 

;.    ■-,  tjr.':  rj  v;:H-'  i-:  '.  r.-r':  .<::.:■  x  •-■  ■•• 
structure  of  size  0(n^  altogether,  while  giving  a  query  time  of  OQogn).  We  can  do  better,  how- 

.•  ''■•  cV'&Ji.  T'l  jfu—^:- •^.  n  '  .? 
ever. 

As  before,  by  forming  a  spanning  tree,  etc.,  we  obtain  a  sequence  of  0{rv')  subsets  of  X  such 
that  adjacent  sets  disagree  on  just  one  point  and  every  solution  to  the  half-space  problem  lies  in 
the  sequence  of  sets.  We  divide  this  sequence  into  (9(n-**)  collections  of  n-"'  adjacent  sets;  for 
each  collection  we  store  a  list  of  the  0{n)  points  contained  in  all  n^"*  sets  in  this  collection.  We 
next  divide  each  collection  into  n*  collections  of  n'~^  adjacent  sets;  for  each  collection  we  store  a 
list  of  the  0(«^~*)  points  contained  in  all  n^"^  sets  in  this  collection,  and  which  have  not  been 
stored  in  some  other  list  associated  with  a  supercollection  of  the  current  collection.  We  then 
divide  each  collection  into  n*  collections  of  n^~^*  adjacent  sets,  and  so  on,  until  we  obtain  collec- 
tions having  exactly  one  set.  For  each  such  collection  we  store  a  list  of  the  remaining  0{n*)  points 
that  have  not  been  stored  in  a  list,  associated  with  a  supercollection  of  the  current  collection.  For 
each  of  the  lists  we  have  aeated  we  build  a  data  structure  to  answer  3-dimensional  half  space 
queries.    The  complexity  of  this  solution  to  the  wedge  problem  is  Oin^**,-^^).    (By  using  the 
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more  complicated  dividing  of  section  4,  we  can  get  0(lognlog(l/€))  time,  etc.,  but  for  simplicity  in 
this  section  we  ignore  this  slight  improvement.) 

We  now  complete  the  solution  of  the  cone  retrieval  problem.  Consider  an  axis  in  the  base 
plane  (the  plane  in  which  the  comer  of  the  cone  lies)  through  the  comer  of  the  cone,  parallel  to 

r  .r  .J  ■'.'       .;■■, 

the  jT-axis.  And  consider  rotating  a  plane  about  this  axis.  We  see  there  are  three  positions  of  this 
plane  for  which  one  of  the  three  lines  defining  the  cone  lies  in  this  plane.  These  three  positions 
divide  the  point  set  into  four  subsets;  we  need  to  perform  a  wedge  query  on  two  of  these  subsets. 
As  each  subset  is  defined  by  at  most  two  planes,  for  each  position  of  the  comer  there  are  O(n-) 
subsets  of  points  that  may  be  searched.  We  show  below  how  to  economise  on  storage;  for  now 
we  continue  with  the  solution. 

The  next  issue  is  to  find  the  number  of  rotational  orders  of  the  n  pomts  with  respect  to  a 

^m:'  -^.^  sc'  c-  ^;jr,ii.p  3-.J-:  - 

plane  rotated  about  an  axis  parallel  to  the  jc-axis^  the  axis  lying  in  the  base  plane.   (And  hence  the 

■  ■:  '.^  ■■'.  ■;    t:nn  •';  jp  e  it 

number  of  times  we  need  to  store  O(n-)  search  structures  for  subsets  of  the  points.)  An  equivalent 
problem  is  obtained  by  projecting  onto  the  >'-r  plane;  it  is  then  dear  there  are  0{rr)  orders.  This 
leads  to  an  0(n^*'"^«)  storage  requirement  for  the  cone  retrieval  problem.  We  show  below  how  to 
reduce  this  somewhat. 

First,  we  consider  the  initial  part  of  the  solution:  rotating  a  plane  about  an  axis  through  a 
specific  comer.  The  n  points  are  encountered  in  a  specific  order.  We  proceed  by  a  technique 
analogous  to  repeatedly  dividing  a  list  into  two  eqtial  sized  sublists.  Then  a  query  for  the  points  in 
a  given  interval  of  the  list  is  reduced  to  the  union  of  O(logn)  queries  on  disjoint  sublists.  Instead 
we  divide  this  order  into  n'  sets  of  n~*  points;  each  of  these  is  in  turn  divided  into  n*  sets  of  n^~^ 
points,  and  so  on.  We  build  a  data  structure  for  carrying  out  a  wedge  query  for  each  of  the  n-* 
possible  sets  of  0(n)  points  created  by  unions  of  adjacent  sets  among  the  sets  of  n^"*  points,  and 
likewise  for  the  sets  obtained  in  the  recursive  division.  Recall  a  cone  query  is  the  union  of  two 
wedge  queries.  Each  wedge  query  becomes  the  union  of  at  most  1  query  on  a  set  of  0{n)  points, 
2  queries  on  sets  of  0{n'~*)  points,  ...,  and  2  queries  on  sets  of  O(n')  points.   So  for  a  given  rota- 
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tional  order  on  n  points  the  query  time  (for  a  cone  query)  is  0(  °^^)  and  the  storage  require- 

raent  is  0(/i^^^').  ■"'"   ■  '"'   ^-  '"     - 

Second,  we  consider  the  latter  part  of  the  solution.  We  note  that  'adjacent'  rotational  orders 
differ  on  exactly  2  points.  So  we  divide  the  rotational  orders  into  0(n^**)  sequences  of  n'~«  adja- 
cent orders;  for  each  sequence  we  store  one  list  of  the  points  on  which  all  the  orders  in  the 
sequence  agree.  For  this  list  we  build  the  data  structure  described  in  the  first  part  of  the  solution. 
Then,  we  recursively  divide  the  sequences  of  size  n^~*  into  n'  sequences  of  size  n^"-',  etc.  .  Now, 
a  search  requires  one  to  examine  the  data  structures  for  lists  of  length  0(n),  0(n^~'),  ...,  0{n*), 

respectively,  taking  time  0{-^^).  The  space  requirement  is  0(n^^**). 

pi;i±t    sa:  vo    r:.  Vc  vl,  .: 

Aligned  Tetrahedral  Retrieval  "^  ^  "  '^^  '^^-^  --^'-^  ='       ^'  '"'^^ 

■-.:■  •Jcr.rrl.'TTCO 'jk: -.?-jO  "'r  i,.^    .r. 
Fu^t,  we  define  a  subproblem:  the  aligned  wedge  query.    This  is  a  wedge  query  with  the 

mtersection  of  the  two  planes  defining  the  half-sp"aces  lying  in  the  base  plane.  A  little  thought 
shows  that  the  solution  for  the  half-space  problem  in  3  dimensions  also  can  be  used  for  this  prob- 
lem; the  complexity  is  0(/)^log^).  '  "'"- -  -^^  "f-  '•^'  •'•''  ■■■■    ~ 

The  aligned  tetrahedral  query  is  the  intersection  of  two  aligned  wedge  queries,  where  the 
aligned  wedges  are  on  separate  parallel  planes  (z=0  and  z=l,  for  example).  It  is  not  hard  to  see 
that  an  aligned  tetrahedral  query  is  defined  by  two  points  in  each  of  the  parallel  planes.  A  simple 
method  would  be  to  keep  every  solution  for  the  first  aligned  wedge  query,  and  for  each  solution 
to  build  a  data  structure  for  carrying  out  the  second  aligned  wedge  query.  The  space  requirement 
would  be  C>(n-°).   We  improve  this  to  a  complexity  of  0(/i'^«,-^^). 

The  idea  is  to  reduce  the  number  of  different  solutions  to  the  first  query  that  we  store,  by 
only  storing  the  common  parts  of  such  solutions  once.  This  reduction  comes  in  two  stages.  The 
first  is  to  reduce  the  number  of  lists  on  which  the  first  wedge  query  may  be  carried  out.  Recall 
that  there  are  0{n*)  different  lists  we  could  search.   We  can  form  a  linear  sequence  of  these  lists. 
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of  length  0(n'*),  possibly  with  some  duplicates,  but  incliidiflg  every  list,  such  that  adjacent  lists 
only  differ  on  two  items.  (Take  spanning  trees,  etc.  .)  We  divide  this  sequence  into  0(n^^')  seg- 
ments of  n'~'  adjacent  lists  each.  For  each  segment  we  extract  a  list  of  the  0(n)  points  on  which 
the  n'~*  lists  in  the  segment  agree.  The  segments  are  divided  into  n*  subsegments  each  holding 
n'~-'  lists.   For  each  subsegment  we  extract  a  list  of  the  0(n^~*)  points  on  which  the  n'~~'  lists  in 

. '.'     :".    ■'!:.:        ...  •■       t 

the  subsegment  agree,  which  have  not  been  stored  in  a  list  associated  with  a  larger  segment  (con- 
taining  the  current  subsegment).  Next  we  divide  each  subsegment  into  n*  subsegments  each  hold- 
ing n^"^*  lists,  etc.  .   For  each  of  the  lists  extracted  in  this  way  we  still  have  to  be  able  to  perform 

'  ,      ■;      :.i  lo  x-'W  ..:.     :z.r^       -?  r. 
the  first  wedge  query. 

:'■*"  'rt)0  ?r  f.f.T.  ••:'pc.  -..■x-  ■  - 
The  second  reduction  comes,  as  already  seen,  by  the  technique  analogous  to  repeatedly 

dividing  a  list  into  two  equal  sized  sublists.    We  divide  each  list  into  n*  equal  sized  sublists.    For 

each  set  of  equal  sized  sublists  we  form  the  n-^  possible  combinations  of  adjacent  sublists.   Then 

any  query  is  a  union  of  0(l/£)  disjoint  lists,  at  most  one  of  length  0(n),  two  of  length  0(n-"'), 

two  of  length  C>(n"~-'),...,  and  two  of  length  0(n*). 

To  each  list  extracted  by  the  first  reduction  we  apply  the  second  reduction.  We  note  that  the 
result  of  the  first  aligned  wedge  query  is  a  union  of  some  of  these  doubly  reduced  lists  (in  fact 
0(1/ 1~)  of  them).  Thus  to  perform  the  second  wedge  query  it  suffices  to  build  a  data  structure  for 
carrying  out  wedge  queries  for  each  of  the  lists  obtained  by  applying  the  second  reduction.  Care- 
ful computation  shows  the  complexity  is  0{n''*^*,  °^"). 

Tetrahedral  Retrieval 

It  remains  to  reduce  the  tetrahedral  query  to  a  union  of  based  cone  queries  and  aligned 
tetrahedral  queries.  We  proceed  in  a  similar  way  to  the  method  for  the  triangle  retrieval.  We  cut 
the  space  by  parallel  planes,  such  that  between  adjacent  planes  there  are  0(n^~*)  points.  There 
are  two  cases:  either  a  subproblem  contains  a  comer  of  the  tetrahedron  or  it  does  not.  In  the 
former  case,  where  we  have  a  query  on  a  point  set  of  size  n"~*,  we  continue  recursively  dividing. 
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The  latter  case  is  the  harder.   There  are  two  possibilities:  either  the  two  planes  defining  the  subset 
of  points  both  intersect  the  tetrahedron  in  a  quadrilateraToflHejrboth-ijiiexiect  the  tetrahedron  in 


a  triangle. 


f^. 


The  latter  case,  both  intersecting  in  triangles,  we  ca^  a  pri^.  4t' is  illustrated  below. 


A  C 

The  prism  query  can  be  seen  to  reduce  to  the  union  of  two  based  cones  ABCF  and  DEFB  and  an 

aligned  tetrahedron  ABDF. 

:■■'   arm  t;  t'jr.r  rt  '.'.-  ':  "  -Bii  ■ 
The  former  case,  both  intersecting  in  quadrilaterals,  is  illustrated  below. 


WLOG  assume  that  ^  <  -ttt.    Let  ^  be  the  point  on  BC  such  that 

flC        FG  OK        t(j 

query  is  the  union  of  the  prism  EFGKAB  and  a  query  consisting  of  two  parallel  faces,  one  a  trian- 
gle and  one  a  quadrilateral.   It  is  illustrated  below. 
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A . 

But  this  is  just  the  union  of  the  three  based  cones  ADHC,  AKCG,  EHGA  together  with  the  aligned 
tetrahedron  A G//C.   We  know  how  to  solve  all  of  these.  ^ 

So  the  solution  to  the  tetrahedral  query  has  complexity  Oitf'^^,^^)^ 


■:■■  \^     .<    -tC-  b^""'  0W3  'L 
7.  Other  Problems 


retrieval:  The  query  is:  'wiiat  points  lie  inside  a  ring  with  inner  radius  r^  and  outer 
radius  r,?'.   This  is  an  extension  of  the  drde  retrieval  problem,  mentioned  above.   We  achieve  a 


complexity  of  C>(n"*,Iogn).  •:  .^:.  -      yi,      ~.     ^ 

We  divide  the  plane  into  regions  using  the  perpendicular  bisectors  of  lines  joining  pairs  of 
points,  forming  a  region  graph  of  size  O  («■*).  At  any  location  within  a  given  region  the  ordering 
of  the  points  by  distance  is  the  same;  we  keep  this  order  as  a  list.  So  we  have  0(n*)  lists  of  length 
n,  ,with  lists  corresponding  to  neighboring  regions  differing  on  exacdy  two  points.  These  lists  do 
not  lie  in  a  sequence;  however,  using  techniques  already  presented,  they  can  be  replaced  by  a 
sequence  of  0{n'^)  lists,  which  includes  every  list  and  has  the  property  that  adjacent  lists  differ  on 
exactly  two  points  (Take  a  spanning  tree,  etc.).  Using  the  method  of  [Co]  we  can  store  the  0(n'^) 
lists  in  space  0(n'^)  while  allowing  a  range  query  on  the  list:  namely,  report  the  points  between  dis- 
tance r,  and  r,  from  the  given  location.  The  range  query  runs  in  time  0(logn+s),  where  s  is  the 
number  of  points  reported.  Fmally,  to  determine  which  list  to  search  we  need  a  point  location 
algorithm,  which  we  apply  to  the  region  graph  of  size  0(n'*).  K  we  use  either  Kirkpatrick's  [K]  or 
Cole's  [Co]  solution  we  achieve  a  complexity  of  0{n*,\ogn). 
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Point  closest  to  a  line:  The  query  is:  Vhidi  point  is  closest  to  a  given  line?'.  We  obtain  a  solution 
of  complexity  0(n-,logn).  We  say  those  lines  whidi  separate  the  points  into  the  same  two  sets 
form  a  channel.  There  are  0(n-)  channels.  For  a  lia&ia:a  givea-diannel,  determining  the  nearest 
point  is  just  a  matter  of  searching  two  convex  hulls  (the  convex  hulls  of  the  two  sets  into  which 
the  points  are  divided  by  the  line).  ActuaHyi,  for^^e^  pair  of  these  convex  hulls  we  only  need  to 
keep  two  segments;  those  determined  by  the.ifflnunonlangeats_ (with  respect  to  die  two  convex 
hulls)  passing  through  the  channel.  On  traversing  such  a  segment,(Mdthjrespect  to  a  line  passing 
through  the  corresponding  channel,  we  find  the  minimum  distance  from  a  point  on  the  segment  to 
the  line  is  unimodular.  Hence  a  Fibonacci  seardi  allows  the  nearest  point  to  be  found  in  time 
O(logn).   Notice  each  line  segment  joii'ing  tw^-jjoints  lies' on  exactly  twq^.convex  hulls;  thus  stor- 


ing the  convex  hulls  takes  space  0{rr). 

The  method  for  determining  the  channel  uses  the  first  and  second  techniques  of  this  paper. 
We  draw  a  baseline;  for  any  point  on  the  baseline  CX")  of  the  diannels  are  Sisible'.  In  fact,  for  a 
fixed  point  on  the  baseline,  the  range  of  a  visible  channel  is  defined  by  two  vertices,  one  on  each 
of  the  convex  hulls  defining  the  chaimel.  We  deduce!  tHaTit  siSfices  to  store  the  data  stnK:ture 
used  for  the  half  plane  retrieval  problem,  if  in  addition,  for  each  range  we  keep  a  pointer  to  the 
appropriate  channel. 

Spherical  retrieval  hi  higher  dimensions:  The  annular  retrieval  can  easily  be  generalized  to  sphere 
retrieval  in  3-dimensions,  and  also  to  higher  dimensions.  The  complexities  are  still  polynomial, 
but  the  polynomial  become  less  modest  as  the  dimension  increases,  mostly  because  of  the  space 
required  for  point  location. 

By  the  remark  that  the  method  of  [W,Y1]  apparently  fails  for  dimensions  >5,  it  follows  that 
their  technique  fails  for  spherical  retrieval  of  dimension  >4. 


Problem 

Previous  Result 

Our  Result 

triangle  retrieval 

;■'   n  .       r 
(XnMcgn)  [EKM] 
CKn,n°^  [W] 

0(ni  lofn) 

annular  retrieval 

not  csmsidered 

O(^^log^) 

drcle  retrieval 

CXn,nO^«)[Yir 

for  fixed  radius 
0(nM6gn)  [BM]->i 

O(nlog3n,logn) 

.probabilistically 
'O(«log2n,logn) 

point  closest  to 
a  line 

not  considered 

CXnMogn) 

half-space  retrieval 
in  3  dimensions 

0(n,n°-^'[Yr] 

-,0  e- Ij^..     . 
CXnMogn) 

half  space  retrieval 
in  d  dimensions 

not  considered 

0(^n^-\2''logn) 

sphere  retrieval 

not  considered 

O(nMogn) 

tetrahedral  retrieval 

not  considered 

0(.-',i^) 
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